Disorder-Induced Resistive Anomaly Near Ferromagnetic Phase Transitions 



Carsten Timm/ M. E. Raikh,^ and Felix von Oppcn 

'institut fiir Theoretische Physik, Freie Universitdt Berlin, Amimallee I4, 14195 Berlin, Germany 
'^Department of Physics, University of Utah, Salt Lake City, UT 84112 
(Dated: August 27, 2004) 

We show that the resistivity p{T) of disordered ferromagnets near, and above, the Curie temper- 
ature Tc generically exhibits a stronger anomaly than the scaling-based Fisher-Langer prediction. 
Treating transport beyond the Boltzmann description, we find that within mean-field theory, dp/dT 
exhibits a, \T — T'c|~^/'^-singularity near Tc- Our results, being solely due to impurities, are relevant 
to ferromagnets with low Tc, such as SrRuOa or diluted magnetic semiconductors, whose mobility 
near Tc is limited by disorder. 
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Introduction. — It was first observed by Gerlach 
that the resistivity of itinerant ferromagnets exhibits an 
anomalous temperature dependence in the vieinity of the 
Curie temperature Tc. This feature was later reproduced 
with much higher experimental accuracy by Craig et al. 
as well as others ■ Dating back to the seminal works 
of de Gennes and Friedel Q as well as Fisher and Langer 
13, this resistive anomaly is conventionally explained in 
terms of coherent scattering of carriers by large blocks 
of spins whose size is determined by the magnetic corre- 
lation length £,{T). As the temperature approaches Tc, 
^(T) diverges, making the scattering more cfhcient. 

de Gennes and Friedel studied the resistive anomaly 
of ferromagnets within mean-field (MF) theory. They ar- 
gue that due to critical slowing down, spin fluctuations 
can be treated as effectively static. The ensuing result 
for the coherent transport scattering rate from s^in fluc- 
tuations above Tc is succinctly summarized by 



— = — I dx 



t + (kpax) 



(1) 



Here, the scattering rate is normalized to the rate 1 /rg for 
incoherent scattering. The denominator of the integrand 
arises from the conventional MF (Ornstein-Zernike) cor- 
relator for spin fluctuations, with t = [T — Tc)/Tc denot- 
ing the reduced temperature, a is a microscopic length 
of the order of the lattice constant. The integration vari- 
able X denotes the transferred momentum in units of the 
Fermi wave vector kp- The numerator incorporates the 
usual factor 1 — cos 9 (with 9 the scattering angle) in the 
transport scattering rate. The integral in Eq. yields 
To/r(i) = To/r(0) -f {l/^{kFaf)t\nt, implying a singu- 
larity of the resistivity of the form p{T) = po — bt\n{l/t) 
with 6 > when approaching Tc from above. 

Fisher and Langer noticed that this singularity 
emerges from the lower limit of the integral in Eq. 
while the body of the integrand is dominated by large 
wave vectors. Within MF theory, the large-wave-vector 
behavior of the spin-spin correlator is non-singular. The 
central assertion of Ref. is that there exists a singular 
contribution to this correlator at large wave vectors when 



going beyond the MF approximation, with the singularity 
governed by anomalous dimensions. Fisher and Langer 
conclude that while below Tc, the predictions of Ref. [J| 
are essentially correct, there is no singularity within MF 
approximation when approaching Tc from above. 

This conclusion rests on their important physical ob- 
servation that Eq. ^ becomes inapplicable for t <C 
(f/a)^, or equivalently, when the correlation length ^(T) 
exceeds the mean free path £ (due to phonon or impu- 
rity scattering). The reason is that the carriers can be 
viewed as plane waves only over distances shorter than £ 
and thus, they are no longer susceptible to the order in 
the spin configuration beyond £. 

Later, the ideas of Refs. 0,01 were extended to include 
realistic features of ferromagnets and to the critical 
behavior of other quantities such as the spin-flip scatter- 
ing rate 0j H, ■ The effect of a flnite mean free path on 
the resistivity was studied by a number of approaches, 
ranging from replacing the 5-function in the golden rule 
by a Lorentzian [Sj to smearing the Ornstein-Zernike cor- 
relator in order to eliminate the pole [i(tI |. 

All of these approaches are based on the Boltzmann- 
equation formalism. The prime message of this paper 
is that in the small-t limit, when the correlation length 
^(T) exceeds the mean free path £, the Boltzmann ap- 
proach fails. The reason for this is that for £,{T) ^ £, 
the smooth variations of the magnetization (on the scale 
^(T)) are "explored" by diffusing carriers. By contrast, 
the Boltzmann approach prescribes to treat scattering 
from both short-range impurities and the smooth varia- 
tions of the magnetization on equal footing, i.e., to add 
their partial scattering rates. 

The consequences of going beyond the Boltzmann ap- 
proach are drastic. In fact, as demonstrated below, in- 
stead of smearing the resistive anomaly, impurities cause 
a much stronger singularity, even within MF theory. 
Quantitatively, we find d^p/dT^ ^ t~^/'^ sufficientlyclose 
to Tc, as opposed to d^p/dT^ - t-^ for ^(T) < ^ g. 

Our reasoning goes as follows. The adequate descrip- 
tion of electric transport on scales larger than the phase- 
breaking length is a network of resistors, made up of 
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FIG. 1: Schematic representation of resistor network describ- 
ing disordered ferromagnets when the correlation length £,{T) 
exceeds the phase-coherence length L^. Each block of size 
constitutes a resistor with coarse-grained magnetization 
(arrows) and random impurities (dots). 

cubes of size L^, as illustrated in Fig. ^ This network 
is inhomogeneous due to the different spin and disorder 
configurations in each cube. Contrary to the Boltzmann 
prescription, it is essential to first compute the effective 
resistivity of the inhomogeneous network and to perform 
the disorder and thermal averages only as the last step of 
the calculation. Then, correlations between distant spins 
affect the resistance through the inhomogeneous current 
and field distribution over the network, leading to a true 
singularity for i ^ even within MF theory. 

To illustrate the importance of performing the impu- 
rity and thermal averages at the last step of the calcu- 
lation, consider a minimal model of two macroscopic re- 
sistors in sequence. These resistors differ in both im- 
purity and spin configurations. Within the Boltzmann 
approach, both resistances are equal to 1/a upon con- 
figurational and thermal averaging, yielding a total re- 
sistance of 2/ct. However, for the actual distribution of 
impurities and spins, their conductivities differ, so that 
(7i = a + 5a /2 while CT2 = a — 5a /2. Then, the effec- 
tive resistance becomes {ai + a2) / aia2 — 2/a + 5a'^/2a^ 
(assuming 5a <^ W). This involves an additional term 
5a^ /2a^ with non-zero average. 

Effective conductivity of an inhomogeneous medium. — 
Remarkably, the effective resistivity p^s = 1/aeS can be 
computed for an arbitrary realization a(T) of the local 
conducti vity , provided that the relative variation in a{r) 
IS weak hj. To see this, we decompose the current, 
the conductivity, and the electric field into averages jo, 
(To, and Eo and spatially fluctuating contributions <5j(r), 
da{r), and (5E(r). From Ohm's law, we have 

jo = rToEo + (<5a(r)JE(r)), (2) 
5j(r) = ,5a(r)Eo + ao5E(r). (3) 



Here, the brackets denote a spatial average. Combining 
the continuity equation V • (5j = and Maxwell's equa- 
tion V X (5E = with Eq. (Q, one obtains (SE(q) = 
— qq- E(5f7(q)/(To, where q denotes the unit vector in the 
direction of the wave vector q. Inserting this into Eq. 
(O, one finds for the effective macroscopic conductivity 
acff (defined by jo = (JcffEo) in three dimensions 

{[6a{r)]^) 

(Jeff = 0-0 t: • (4) 

3(70 

Remarkably, this result is independent of the geometry of 
the conductivity variations. To illustrate this, consider 
a sample with a 50-50 random mixture of domains of 
conductivities ai and a2, where \ai—a2\'t^ai,a2. In this 
case, Eq. Q yields acS = ao — {ai — a2)^/12ao , which is 
completely independent of the arrangement of domains. 
It is interesting to remark that in two dimensions, one 
can find an exact, geometry- independent expression for 
such two-phase systems which is valid for arbitrarily large 
inhomogeneities \a2 — ai\/{ai + (T2) E3- 

Effective conductivity due to spin fluctuations. — The 
conductance of a phase-coherent sample exhibits random, 
sample-specific, and reproducible variations as a func- 
tion of external parameters such as the Fermi energy Ep- 
These conductance fluctuations arise due to interference 
between different elastic scattering paths of a carrier that 
diffuses through the sample Thus, the conductance 
g{r, Ep) varies from block to block because of their dif- 
ferent impurity configurations. The fluctuations in the 
conductivity entering Eq. can then be expressed as 

Saiv)^5g{r,EF;mir))/L^. (5) 

We assume that the system is so large that domains with 
any magnetization combined with any disorder realiza- 
tion appear. Using this ergodicity assumption we replace 
the spatial average in Eq. by independent disorder 
and thermal (magnetization) averages. 

Two spin subbands. — In our case, the external para- 
meter is a vector, namely the magnetization ni(r). To 
proceed, we first consider the simplest case in which the 
dominant effect of the impurity spins on the carriers is 
an effective Zeeman field arising from the exchange inter- 
action, which is proportional to the magnetization m{r). 
Then the coarse-grained magnetization m(r), which is 
averaged over each cube of size L^, can be incorporated 
via equal, but opposite energy shifts ±am{r) for spin-up 
and spin-down carriers, and 5a becomes a sum of contri- 
butions from the two spin projections, 

5a{r) = [5g^ir,EF) + 5g^{r,Ep)]/L^ 

= [6gir,E+{r)) + 5gir,Epir))]/L^. (6) 

Here, g{r, Ep) is the spinless conductance of the cube 
at position r and we have introduced the shifted Fermi 
energies Ep ~ Ep±am, see Fig.|5| We used that carriers 
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FIG. 2: Sample-specific variation of g{Ep) for a block of the 
resistor network (UCF) in the absence of spin. The conduc- 
tances for each spin direction are obtained by including equal, 
but opposite exchange- induced Zeeman shifts of the Fermi en- 
ergy. As indicated, this leads to a difference in the conduc- 
tances for the two spin directions. 



with both spin projections are scattered from the same 
impurities. Substituting Eq. © into Eq. we obtain 



Pcff -P0 = -fl {([<5,9(r,i?+(r))]2) + ([<5g(r, i?-(r)r) 
+2 {dg{r, Ep + am(r)) Sg{r, Ep - am{r)))} . (7) 



The first two terms on the rhs describe universal conduc- 
tance fluctuations, which are independent of the Fermi 
energy and thus of m(r) 0|. Therefore, the only t- 
dependence comes from the conductance correlator. 

It is well known that this conductance correlator 
is a function F{x) of the dimensionless ratio x = am/ Ec, 
where Ec = D/L'^ is the inverse diffusion time through 
a block of size (D denotes the diffusion constant). 
In three dimensions, the asymptotic behaviors of F{x) 
are given by [l^ F{x) = F(0)(1 - Cix"^) for a; < 1 and 
F{x) = F(0)C2a;"^/^ for x > 1, where Ci and C2 are 
constants of order unity. The remaining step is to substi- 
tute F{x) into Eq. jT)) and to perform the thermodynamic 
average over m, using the MF distribution for t > 0, 

7'[m(r)] cx exp J cPr [a^ Vjin • Vjin-f i m^] 

(8) 

Here, ca^ is a spin stiffness. Using this distribution, one 
finds for the magnetization fluctuations 



(m2(r)) 



(9) 



The restriction in the range of the q integration accounts 
for the fact that we are computing fluctuations of the 
coarse-grained magnetization. The appearance of t in 
the denominator of Eq. Q manifests the fact that the 
net strength of fluctuations grows when approaching 
Tc- This is in contrast to the conventional origin of the 
t-dependence, namely the spin correlations. Perform- 
ing the integration in Eq. ©, we obtain (m^(r)) ~ 



(3T,/2^2^a2) - (l/e(T)) arctan(e(T)/L^)] in 

terms of the MF correlation length ^(T) ~ a/^/t. 

As the next step, we replace TO(r) in Eq. O by 
(m^(r))^/^. Using the distribution ©, this procedure 
can be shown to be exact in the limits of small and 
large x. Expanding the result in the small parameter 
Lc/,/£,{T) = L^^/t/a^ we readily obtain 



Peff — PO 



F{xo)- F'{xo)xo 



4a 



(10) 



where xq = {alEc){iTj2iT'^ca?L^f/'^. The t-dependent 
part of pcs is given by the second term on the rhs 0| . 

So far, our model completely disregards spin-orbit 
(SO) coupling which is present in the vast majority of 
ferromagnets. Below, we incorporate SO coupling into 
the calculation of the effective resistivity. 

Spin-orbit coupling. — In the presence of SO coupling, 
the variance {[6g{r, Ep;m.{r))]'^)i„^p increases by a fac- 
tor of two when applying a sufhcicntly strong Zee- 
man field In terms of the relevant dimension- 
less measure of the exchange-induced Zeeman field y = 
amTso (here Tgo is the SO time), we can then write 
USQir,EF;mir))]^)imp = Hiy) with H{oo)/H{0) = 2 
[1^. Combining Eqs. lO and substituting H{y), and 
proceeding as in the derivation of Eq. Hl()|l . we obtain 



PcS - Po 



3Ll 



Hiyo)^H'iyo)ya 



4a 



(11) 



where yo — aTso{3Tc / 27r^ ca^ L fj,)^ . Eq. ((TT|l assumes 
that the SO length iso = (-Dtso)^^^ is smaller than L^. 

Discussion. — By going beyond the Boltzmann descrip- 
tion in calculating the critical behavior of the resistivity, 
we have expressed poff through mesoscopic characteris- 
tics. Recall that our results, Eqs. (|10|l and were 
obtained within MF theory so that the singularity in the 
t-dependence of pQ on the Ihs of these equations, aris- 
ing within the Boltzmann formalism, is suppressed by 
impurities, since ^{T) > £. Eqs. (^J and ((TT|l show 
that in addition to this suppression, impurity scatter- 
ing leads to a singularity which is much stronger 
than the de Gennes-Friedel result t\n{l/t) and the large- 
wave-vector Fisher-Langer contribution. This singular- 
ity, which in essence is governed by Kirchhoff 's laws, con- 
stitutes our central result. To resolve the anomaly on top 
of a monotonous phonon contribution to p, it is custom- 
ary to consider (Pp/dT^. Then, our disorder-induced MF 
anomaly becomes cPp/dT^ ~ t^^^^ . In a log-log plot, the 
slope is —1 in the Boltzmann regime £,{T) ^ £ and —3/2 
for the disorder-induced anomaly, as shown in Fig. |31 

According to Eqs. (|10|l and Hll() . the sign of the 
anomaly is governed by the sign of either F'{xo) or 
H'{yo)j depending on the strength of SO coupling. Since 
F'{xo) < while i?'(yo) > 0, the disorder-induced 
anomaly corresponds to a decrease of d^p/dT^ when ap- 
proaching Tc from above in the absence of SO coupling 



4 




,' without SO coupling 



FIG. 3: Resistive anomaly within MF theory for ^ I 
(schematic) with (full line) and without (dashed line) SO cou- 
pling. The anomaly is described by the de Gennes-Friedel 
mechanism for ^(T) ^ ^, while the disorder-induced mecha- 
nism of this paper dominates closer to Tc where ^(T) ^ L^. 
When ^ I, there is an additional intermediate regime. 
Inset: anomaly with SO coupling in a log-log plot. 

and to an increase in the realistic case of strong SO cou- 
pling. The difference in signs comes about because the 
Zeeman field suppresses the correlator in Eq. (O, while 
it increases the UCF for strong SO interactions. 

The magnitude of the disorder-induced anomaly is con- 
trolled by F(0) and H(0), which are the variances of the 
conductance, (((5g)^)imp, cf. Eqs. ifTUI) and lfTT)l . This 
quantity assumes different values in different regimes 
which are defined by the relations between the rele- 
vant lengths, namely L^, the thermal length = 
{D/Tf/^, and the spin-flip length £s- In the simplest 
case, when is the smallest length, LT,is, we 

have {{Sg)'^)inip ~ (e'^/h)'^ 0. If is larger than Lt 
or £s, then ((5(7)^)imp is suppressed below the universal 
limit, {iSg)^)imp ^ (eV/i)2[min{LT,4}/i0] HI. 

The disorder-induced anomaly proposed in this paper 
is most relevant to ferromagnets with high resistance and 
low Tc since in such systems (i) the mean free path is 
dominated by impurity scattering [l^ and (ii) the phase- 
coherence length exceeds the mean free path at Tc- Nat- 
ural candidates for disordered low- Tc ferromagnets are 
SrRuOs 01 which belongs to the class of poor metals 
1^ as well as diluted magnetic semiconductors (DMS) 
13. l2oj which lately attracted considerable attention in 
view of possible spintronics applications Indeed, 
both types of materials show resistive anomalies which 
differ significantly from the predictions of Fishcr-Langer 
theory. Metallic samples of DMS exhibit pron ounced 
maxima near even in p vs. T HHllillillillilli^, 
which are unrelated to the T = metal-insulator transi- 
tion imi^. 

In closing, it is interesting to point out that our prin- 
cipal result, namely the enhancement of the resistive 
anomaly by disorder, can be viewed in perspective of the 



enhanced coupling of the spin fluctuations to the carriers, 
expected due to the diffusive carrier dynamics. 
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